HEPHY-PUB 693/98 
UWThPh-1998-38 
|hep-ph/990439T 
April 1999 



QUALITY OF VARIATIONAL 
TRIAL STATES 



Wolfgang LUCHAQ 

Institut fiir Hochenergiephysik, 
Osterreichische Akademie der Wissenschaften, 

Austria 



Franz F. SCHOBERLfi 

Institut fiir Theoretische Physik, 
Universitat Wien, Austria 



Abstract 

Various measures for the accuracy of approximate eigenstates of semibounded self-adjoint 
operators H in quantum theory, derived, e.g., by some variational technique, are scrutinized. 
In particular, the matrix elements of the commutator of the operator H and (suitably chosen) 
different operators with respect to degenerate approximate eigenstates of H obtained by the 
variational methods are proposed as new criteria for the accuracy of variational eigenstates. 
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1 Motivation 

A central element in quantum theory is the solution of eigenvalue problems. However, 
usually there is no suitable exact solution to which perturbation theory can be applied. 

A very efficient way to locate the discrete spectrum of some self-adjoint operator H 
bounded from below is provided by the famous Rayleigh-Ritz variational technique 0: 
If the eigenvalues Ek, k = 0, 1, . . ., of H are ordered according to Eq < E\ < E 2 < . . ., 
the first d of them are bounded from above by the d eigenvalues Ek, k — 0, 1, . . . , d— 1, 
(ordered by E < Ex < . . . < E<i-i) of that operator which is obtained by restricting H 
to some d-dimensional subspace of the domain of H, i.e., Ek < Ek, h — 0, 1, . . . , d — 1. 
However, frequently it is not straightforward to quantify how close approximate and 
exact eigenstates are. Thus, we embark upon a systematic study of the accuracy of the 
variationally determined eigenstates of H and suitable measures to judge their quality. 



2 Measures of the Quality of Trial States 

Consider some self-adjoint operator H, W = H, assumed to be bounded from below. 
Suppressing, for the moment, the index k — 0, 1, 2, . . . , let the eigenvalue equation for 

H\x)=E\ X ), (1) 

be solved by some (generic) eigenvector \x) corresponding to some (real) eigenvalue E. 
The Rayleigh-Ritz variational technique yields an upper bound E on this eigenvalue E 
as well as, by diagonalization of the relevant characteristic equation, the corresponding 
vector \tp) in the d- dimensional trial space. There exist several (potentially meaningful) 
measures of the quality of this trial state \ip) which immediately come to one's mind: 

1. The trial state \<p) is supposed to represent — to a certain degree of accuracy — the 
approximate solution of the eigenvalue problem defined in Eq. (TJ). Consequently, 
a first indicator for the resemblance of \<p) with the exact eigenstate \x) would be 
the distance between the expectation value of the operator H with respect to the 
trial state \ip), i.e., between the obtained upper bound E = ((p\H\(p)/{ip\ip), and 
the exact eigenvalue E. However, the precise location of the exact eigenvalue E is 
usually not known. 

2. The natural measure for the resemblance of the Hilbert-space vectors \(p) and \x) 
under consideration is the overlap 

s _ (v\x) , 2 v 
Wi (x\x) 



of the trial state \tp) with the eigenstate \x)- 

Consider the commutator [G, H] of the operator H under consideration with any 
other operator G, where the domain of G is assumed to contain the domain of H. 
Then the expectation value of this commutator with respect to a given eigenstate 
\x) ot H vanishes: 

(X\[G,H]\ X ) = 0. (3) 

Hence, choosing different operators G generates a whole class of operators [G, H] 
each of which may serve to test the quality of a given trial state \(p) by evaluating 
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how close the expectation value (<p\ [G, H] \<p) with respect to \ f) comes to zero. 
This expectation value vanishes, of course, also if, by accident, the state \<p) is an 
eigenstate of G. However, for a given operator G, after having determined \ip), it 
is straightforward to check for this circumstance, for instance, by inspecting the 
variance of G with respect to \(p); the latter vanishes if \ip) is an eigenstate of G. 
Moreover, it goes without saying that an expectation value (cp\ [G, H]\<p) vanishes 
also if the state \(p) is an eigenstate of the commutator [G, H] with eigenvalue 0, 
or even if the state defined by [G, H]\(p) proves to be orthogonal to the state \(p). 

For any self-adjoint operator G, i.e., G^ = G, this commutator is anti-Hermitean, 
which clearly suggests to define a self-adjoint operator C = C' (on the domain of 
H) by [G, H] =: i C. If, for example, G is chosen to be the generator of dilations, 

G = - (x ■ p + p ■ x) , (4) 

the relation (§) is precisely the "master virial theorem" introduced in Ref. || for 
a systematic study of (relativistic) virial theorems 0. In this case, for operators 
H of the form of some typical Hamiltonian consisting of a momentum-dependent 
kinetic-energy operator, T(p), and a coordinate-dependent interaction-potential 
operator, V(x), that is, H = T(p) + V(x), the operator C becomes the "virial 
operator" 

The point spectrum (i.e., the set of all eigenvalues) of the dilation generator (|J) is 
empty; in other words, the dilation generator has no Hilbert-space eigenvectors. 



3 Spinless Salpeter Equation 

Let us apply the above general considerations to the prototype of all (semi-) relativistic 
bound-state equations, the "spinless Salpeter equation," defined by the Hamiltonian 
(in one-particle form, encompassing also the equal-mass two-particle case ||, ||, [6|, [/J) 

H = T + V ; (6) 

here T is the relativistic kinetic energy of some particle of mass m and momentum p, 

T = T(p) = ^Jp 2 + m 2 , 

and V = V(x) is an arbitrary, coordinate-dependent, static interaction potential. The 
spinless Salpeter equation is then just the eigenvalue equation for H, H\xk) = Ek\xk), 
k — 0, 1, 2, . . . , for the set of eigenvectors \xk) corresponding to energy eigenvalues Ef.. 
Analytic upper bounds Ek on these eigenvalues have been given S S S B S H9- 
For the sake of comparison, we focus our interest to central potentials V(x) = V(r), 
r = |x|. Furthermore, in order to facilitate the numerical treatment of the problem, we 
consider the harmonic-oscillator potential 

V(r) = ar 2 , a>0. (7) 

The reason for this particular choice is the following: In momentum space, the operator 
r 2 is represented by the Laplacian with respect to the momentum p, r 2 — > — A p , while 
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the kinetic energy T, nonlocal in configuration space, is represented by a multiplication 
operator. Consequently, exactly for a harmonic-oscillator potential the semirelativistic 
Hamiltonian H in its momentum-space representation is equivalent to a nonrelativistic 
Hamiltonian with some (effective) interaction potential reminiscent of the square root: 

H = -a A p + yp 2 + m 2 . (8) 

The solutions of the corresponding eigenvalue equation may then be found with one of 
the numerous procedures designed for the treatment of the nonrelativistic Schrodinger 
equation. 

For the harmonic-oscillator potential, it is comparatively easy to get a first idea of 
the approximate location of the energy levels Ef. by entirely analytical considerations: 

• On the one hand, every eigenvalue Ek is bounded from above by the eigenvalue 
E k nr of the nonrelativistic counterpart of H: E k < -E*:,nr- 

• On the other hand, every eigenvalue Ef. is bounded from below by the eigenvalue 
E k (m = 0) of the Hamiltonian H corresponding to a vanishing particle mass m: 
E k > E k (m = 0). 



4 The "Laguerre" Trial Space 

As far as the achieved accuracy of the solutions obtained is concerned, the most crucial 
step in all variational games of the Rayleigh-Ritz kind is, for a given operator H under 
consideration, a reasonable definition of the adopted trial subspace of the domain of H . 

For spherically symmetric potentials V(r), a very popular choice for the basis states 
which span the trial space required for the application of the variational technique are 
"Laguerre" trial states, defined in configuration-space representation by JD], [| ||, [7| 



where L k (x) denote the generalized Laguerre polynomials (for the parameter 7) \12\ 



\ T(2£+2f3 + k + l) 



r^eM-^r)L^ p >{2nr)yU^) , (9) 



and D^ m (fi) are the spherical harmonics for angular momentum I and its projection m. 
The trial functions @ involve two variational parameters, fi (with dimension of mass) 
and j3 (dimensionless), which, by the requirement of normalizability of these functions, 
are subject to the constraints \i > and 2/3 > —1. 

One of the advantages of the trial function (|9p is the easy availability of an analytic 
expression for the corresponding momentum-space representation of these trial states. 

For the present investigation, we too employ the "Laguerre" trial states defined by 
Eq. with, for both definiteness and ease of calculation, the variational parameters 
/1 and (3 kept fixed to the values \i = m and (3=1. 



5 Rates of Convergence of the Quality Measures 

Now, let us observe our variational eigenstates, \(p), approaching the exact eigenstates, 
\x), for increasing dimension d of the employed trial space, by comparing the behaviour 
of the various measures for the accuracy of approximate eigenstates discussed in Sec. [| 
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Without doubt, the only genuine "point of reference" of any variational solution to 
an eigenvalue problem is the corresponding exact solution. The exact solution sought is 
computed here with the help of the numerical integration procedure developed for the 
solution of the nonrelativistic Schrodinger equation in Ref. [ |i~3|j . 

Table |I] confronts, for the ground state and the lowest radial and orbital excitations, 
the approximate solutions as calculated with the help of the Rayleigh-Ritz variational 
technique for "Laguerre" trial subspaces of the domain of H of increasing dimension d 
with the exact solutions of the eigenvalue problem for the semirelativistic Hamiltonian 
(^|) with a central interaction potential of the harmonic-oscillator form (0). First of all, 
as discussed in Sec. [| the exact position of any eigenvalue E of our Hamiltonian H is 
confined to a range defined by the nonrelativistic upper bound £"nr and the zero-mass 
lower bound E(m = 0) on this energy eigenvalue E. There are several quantities which 
may participate in a competition for "the best or most reasonable measure of quality:" 

1. The relative error e = (E — E)/E of every upper bound E on the exact energy 
eigenvalue E is, by definition, always nonnegative, i.e., e > 0. 

2. The deviation from unity, a, of the modulus squared of the overlap S of exact and 
variational eigenstates defined in Eq. ([|), a = 1 — \ S\ 2 , is clearly confined to the 
range < a < 1. 

3. The use of the expectation values of the commutators [G, H] with respect to the 
variational eigenstates \ip) is illustrated for the particular example of the dilation 
generator G defined in Eq. (f|), by considering (suitably normalized) expectation 
values (ip\C\(p) of the virial operator C given in Eq. ([5]): 

Mglg) (plp-jr(p)lp) 
{<p\x-£v{x)\<p) {<p\x-£y{x)\<p) 

4. Finally, the normalized maximum difference of the normalized momentum-space 
representations y5(p) and x(p) °f variational eigenstate \ip) and exact eigenstate 
\x), respectively, i.e., the maximum pointwise relative error in momentum space, 
uj = max p [y5(p) — x(p)]/ nrax p x(p) is listed. 

Note that the only measure for the accuracy of approximate eigenstates \ip) which does 
not require any information other than the one provided by the variational technique 
is v, the (normalized) expectation values of the commutator [G, H] with respect to \<p). 
Inspection of Table [I] reveals that v represents indeed a sensitive measure of quality: for 
increasing trial-space dimension d it converges to zero at roughly the same rate as both 
energy and overlap error, e and a, but makes more sense than a pointwise error like u. 

6 Summary and Conclusions 

Various measures for the accuracy of approximate eigenstates of arbitrary (self-adjoint, 
semibounded) operators H have been studied. The vanishing of the expectation values 
of the commutator of H and any other well-defined operator, taken with respect to the 
approximate eigenstates, provides a useful set of criteria for estimating the significance 
of the variational solution. This has been illustrated by considering the commutator of 
the Hamiltonian of the spinless Salpeter equation and the generator of dilations. 
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Table 1 : Characterization of the quality of the variational solution of the eigenvalue problem 
of the semirelativistic Hamiltonian H = \/p 2 + m? + V(r) with harmonic-oscillator potential 
V(r) = ar 2 , for states of radial quantum number n r = 0, 1, 2 and orbital angular momentum 
£ = 0, 1, 2 (called IS, 2S, 3S, IP, and ID in usual spectroscopic notation), obtained with the 
help of our "Laguerre" trial states spanning trial spaces of increasing dimension d = 1, 2, 25, 
by: the nonrelativistic upper bound .Enr and zero-mass lower bound E{m = 0) on the energy, 
the (numerically computed) "exact" energy E, the variational upper bound E on this energy, 
the relative error e of the upper bound, the deviation from unity, a, of the overlap squared of 
exact and variational eigenstates, the (appropriately normalized) expectation values v of the 
virial operator C, and the (normalized) maximum local difference u> of the momentum-space 
representations of exact and variational eigenstates. The physical parameters are fixed to the 
values m = 2 GeV for the particle mass and a = 2 GeV 3 for the harmonic-oscillator coupling. 
A simple entry "0" indicates that the numerical value is closer to than the rounding error. 



Quantity 


d 






State 










IS 


2S 


3S 


IP 


ID 


n r 







1 


2 








I 













1 


2 


£nr [GeV] 




4.12132 


6.94975 


9.77817 


5.53553 


6.94975 


E(m = 0) [GeV] 




2.94583 


5.15049 


6.95547 


4.23492 


5.35234 


E [GeV] 




3.82493 


5.79102 


7.48208 


4.90145 


5.89675 


E [GeV] 


1 


4.21624 






6.50936 


9.77866 




2 


3.92759 


8.10850 




5.24154 


7.18242 




25 


3.82494 


5.79114 


7.48290 


4.90149 


5.89681 


e 


1 


0.1023 






0.3280 


0.6583 




2 


0.0268 


0.4002 




0.0694 


0.2180 




25 








0.0001 








a 


1 


0.09618 






0.36144 


0.65587 




2 


0.02375 


0.43693 




0.09001 


0.34398 




25 








0.00008 








V 


1 


-0.6120 






-0.8328 


-0.9074 




2 


+0.0308 


-0.8666 




-0.5103 


-0.7483 




25 





-0.0001 


+0.0001 








to 


1 


+0.9277 






+0.7541 


+1.0578 




2 


-0.00754 


+2.4577 




+0.3598 


+0.7262 




25 


+0.00003 


-0.0017 


+0.0002 


+0.0004 


+0.0003 



